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Abstract 

The first observation of open b production in ep collisions is reported. An event sample 
containing muons and jets has been selected which is enriched in semileptonic b quark 
decays. The visible cross section a{ep bbX P^X') for < 1 GeV^, 0.1 < 
y < 0.8 is measured to be 0.176 ± 0.016 (stat.) lo'oiy (syst.) nb for the muons to be 
detected in the range 35° < 6^ < 130° and > 2.0 GeV in the laboratory frame. The 
expected visible cross section based on a NLO QCD calculation is 0.104 ± 0.017 nb. The 
cross sections for electroproduction with < 1 GeV^ and photoproduction are derived 
from the data and found to be a{ep ebbX) = 7.1 ± 0.6(stat.) l| |(s?/st.) nb and 
C7(7P ^ bbX) = 111 ± 10 (stat.) (syst.) nb at an average (W^p) ~ 180 GeV, 
respectively. 
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Introduction 



The study of heavy quark production in electron-proton scattering provides an important testing 
ground for QCD. Measurements of open charm production [JI]] at the electron-proton collider, 
HERA, have been shown to be reasonably well described by next to leading order (NLO) QCD 
calculations [0, |], 0, |5|] based on the photon-gluon fusion mechanism. It is essential to extend 
these observations and tests of QCD to the production of heavier quarks. 

In the past, the leptoproduction of open b quarks has been too difficult to observe due to 
the small production cross section. At fixed target experiments only upper limits on the lepto- 
production cross section have been published [§, 0, |Sp []. In this paper, the first observation of 
b quark leptoproduction is presented in the range less than 1 GeV^, where is the four- 
momentum transfered from the electron to the proton. This corresponds to the exchange of an 
almost real photon (photoproduction) i.e. Q'^ ~ 0. 

The data were collected at the HERA collider with the HI detector during 1996 and cor- 
respond to an integrated luminosity of £ = 6.6 pb^^ The production of b quarks is observed 
through their semileptonic decay to muons. The production cross section is measured in a 
kinematic range determined by the properties of the muon from the heavy quark decay and is 
compared with the prediction of a NLO QCD calculation for partons [Q] to which a model for 
the hadronisation and semileptonic decays is added. 



Detector Description 

Only a short description of the detector components which play a major role in this analysis 
is given here; a detailed description of the detector and its trigger capabilities can be found 
elsewhere [Q. Charged particles are measured by two cylindrical jet drift chambers (CJC) 
[|T0|], mounted concentrically around the beam line inside a homogeneous magnetic field of 
1.15 Tesla, yielding measurements of particle momenta in the polar angular rangeQ of 20° to 
160°. Two double layers of cylindrical multiwire proportional chambers [ [TT| ] with pad read- 
out for triggering purposes are positioned inside and in between the two drift chambers. The 



tracking detector is surrounded by a fine grained liquid argon calorimeter [|T2p, covering a polar 
angle range of 4° to 154°, consisting of an electromagnetic section with lead absorbers and a 
hadronic section with steel absorbers. Muons are measured as tracks in the CJC and identified 
by demanding a good link to a track in the instrumented iron return yoke which surrounds the 
superconducting coil. For events from deep-inelastic scattering (DIS) at higher the scattered 
electron is detected by a "spaghetti-type" lead-scintillator-fibre calorimeter (Spacal) [ [131 ] which 
covers the backward region of the detector up to a polar angle of 177.8°. 



Analysis Procedure 

The production of b quarks is detected by looking for muons from the semileptonic decay of 
a b hadron inside jets from the hadronisation of the b quark. A sample of events enriched in 
semileptonic decays of b quarks is selected by demanding that each event has at least 2 jets and 

'EMC also derived a cross section based on the observation of one clear event [|]|. 

1 uses a right-handed coordinate system with the z-axis pointing in the direction of the proton beam (forward) 
and the a;-axis pointing towards the centre of the storage ring. 



that at least one of the jets should contain a muon. The jets are built by using both the calorimet- 
ric and track information and are identified by using a cone algorithm [0] with a cone radius 
r = a/ (Ar])^ + (A0)2 < 1, where is the azimuthal angle and t] is the pseudorapidity and 
with transverse energy to the beam direction, Et, of more than 6 GeV. The muon is required to 
be in the central region of the detector (35° < 6^^' < 130°) and to have a transverse momentum 
with respect to the beam direction, p^, of more than 2 GeV. The selected events have a mean 
track multiplicity of 16.5 and a minimum of at least 5 tracks per event is required. The events 
are triggered by requiring the existence of a muon candidate in coincidence with tracks. The 
trigger efficiency to select these events is determined directly from the data to be 85.5 ± 1.6% 
and is found to be constant for > 2 GeV. To suppress contributions from DIS at higher 
the events which contained an electron candidate with more than 4 GeV energy deposited in the 
Spacal are rejected. To avoid regions of high radiative corrections and to eliminate the remain- 
ingDIS events, the value of y was limited to the range 0.1 <y < 0.8. Here y is g- P/£- P where 
q, P and d. are the four momenta of the virtual photon, incident proton and lepton, respectively. 
It is calculated using the Jacquet-Blondel method JT^]. These conditions limit the data sample 
to photoproduction events with < 1 GeV^; the mean is 6 ■ 10"^ GeV^. In total 927 events 
fulfil these selection criteria. 

Beauty is separated from charm and other backgrounds on a statistical basis using the trans- 
verse momentum of the muon, p^^gp measured relative to the thrust axis of the jet which con- 
tains it. The thrust axis is defined as the direction such that 

T = max 

where the sums run over all the neutral and charged particles in the jet with the exception of the 
muon. Here, pi is the momentum and its longitudinal component with respect to the thrust 
axis. Since the h quark is heavy, it is expected that the spectrum is harder for events orig- 
inating from h quark decays than for decays from lighter quarks. For the determination of the 
beauty cross section, the p'^^.^i spectrum is used for a combined fit of the different contributions 
to the selected data sample. The shapes of the h and c contributions to the are taken from 
the Aroma Monte Carlo [ [T^ event generator, which is based on leading order (LO) QCD ma- 
trix elements, with the photon directly interacting with the partons in the proton and additional 
initial and final state parton showers. 

The method relies on a good understanding of the energy flow within jets, for jets originat- 
ing from heavy quarks, and on a correct determination of the muon misidentification probabil- 
ity. To check that the Monte Carlo simulation describes the energy distribution within jets, 
a data sample enriched in charm quarks is studied, where events with at least one jet and 
at least one D* candidate are selected. The D*^ candidate is identified through the decay 
D*^ D^n^ K^n^n^ and its charge conjugate. This gives a sample of events with 20% 
background which has been verified to have a similar distribution to charm. In each of these 
events jets are identified using the same selection criteria as in the muon sample. The energy 
flow within the jet has been compared to the energy flow as given by the Aroma Monte Carlo 
event generator for charm events with a D* . There is a good agreement between the data and 
the prediction from the Monte Carlo (see figure [1]). From this and many other such comparisons 
it is concluded that the Monte Carlo simulation reliably models the data distributions. 

The selected sample contains muons originating from the semileptonic decays of h and 
c quarks as well as from light hadrons which are misidentified as muons, due to decay or 
hadronic energy leakage (punchthrough). To obtain the background contribution originating 
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Figure 1 : Fractional energy flow within jets as a function of the distance to the jet axis. This fraction is 
the ratio of the energy in an annulus at radius r to the total energy in the jet. The data points (filled dots) 
correspond to a data sample of mainly charm decays and the full histogram to a Monte Carlo simulation 
where only charm quarks are generated. Both the data and Monte Carlo are normalized so that the 
integral of the distributions equals one. 



from misidentifying a light hadron as a muon, it is essential to know the probability V^{p, 0) 
for this to occur. Large samples of single pion, kaon and proton tracks are passed through 
the full detector simulation and the fraction measured as muons is determined. From these, 
the misidentification probability functions Vf^{p, 9), where h = tc, K,p, are parameterised as 
a function of the momentum and the polar angle of the hadron. These functions vary with the 
polar angle but do not exceed 6 • 10^'^ in the case of pions and 2 ■ 10~^ in the case of kaons. For 
protons it is found to be below 2 ■ 10^'^. These probability functions are verified in the data by 
studying Kg and decays, as a source of real pions and kaons, respectively. The distributions 
of fake muons agree both in shape and absolute magnitude (figure ^ with those predicted by ap- 
plying the probability functions to the pions from Kg decays. For kaons from decays, 8 fake 
muons are observed compared to 7.8 predicted. These observations show that the probability 
functions are well determined. 
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Pti[GeV] y^,[deg 
Momentum (a) and polar angle (b) distributions for fake muons in a pion sample from the 
—>■ vr+TT^. The filled points show the fake muon yield as measured in the data. The solid 
histogram gives the estimate of the muon yield as obtained by assigning V'^ip, 9) as a weight to every 
pion and summing the weights over the entire pion sample. The measured yield is 64 fake muons; the 
estimate amounts to 69.8. 
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The knowledge of the 9) functions allows the background from hadrons, which are 

falsely identified as muons in the pi^ spectrum, to be calculated in shape and absolute magni- 
tude from the data. For this purpose, events with at least two jets with Et > 6 GeV each, but 
without the muon requirement, are selected. The expected background in the pt^ ^.^^ distribution 
is calculated from all hadrons which pass the and polar angle requirements. The assign- 
ment of a hadron as a pion, kaon or proton is made using the Jetset[[T7|] specification of the 
fractions of these particles which is in agreement with measurements by the SLD and DELPHI 
collaborations [jTSp. The Jetset pion and kaon fractions are varied by ±0.05, covering the 
largest deviation of Jetset from these measurements. This variation corresponds to about 15% 
relative variation of the kaon fraction. The result of this variation is included in the systematic 
error for this background determination. 

Results 

The observed .^^^ distribution is shown in figure |] together with the b signal and the different 




1 2 3 4 5 



P^.rel [GeV] 

Figure 3: The measured distribution in the data and the fitted sum (solid line) of the contributions 
of beauty (dashed line), charm (dotted line) and the fixed fake muon background (dashed-dotted line). 

background contributions. The signal and backgrounds are obtained from the combined fit of 
the b and c contributions, using the shapes from the Aroma Monte Carlo, together with the 
fake muon background determined as described above. The relative composition of the data 
sample amounts to ft = 50.8 ± 4.9% (beauty), = 22.4 ± 5.0% (charm), and fjake = 25.9% 
(background, fixed). The uncertainty of the latter is included in the systematic error which is 
deduced by changing the pion and kaon fractions as described above. 

The visible electroproduction cross section for b quarks, determined from the number of 
muon events attributed to b quark decays, is measured to be: 

a^isiep bbX fiX') = 0.176 ± 0.016 (stat.) tofit (syst.) nb 



in the kinematic range < 1 GeV^ 0.1 < y < 0.8, p''^ > 2.0 GeV and 35° < < 130°, 
where the first uncertainty is statistical and the second systematic. 

The acceptance is determined from the AROMA Monte Carlo simulation to be 24%, of which 
the efficiency for muon identification is 56% and the probability to reconstruct two jets each 



with Et > Q GeV is 42.9%. The full analysis has been repeated using the Herwig 5.9 [ [T9| ] 
and Rapgap [ PD[ ] Monte Carlo simulations assuming both direct and resolved production of 
the c and b quarks. In addition, different cuts for the jets and different values of the Peterson 
fragmentation parameter, e, for both the c and b fragmentation have been used. The variation 
of the cross section with the different simulations amounts to ±7.1% and this is taken to be 
the systematic error due to the uncertainties in the Monte Carlo. The systematic uncertainty 
due to the muon reconstruction is +6%. The systematic error due to the uncertainty (±4%) 
in the hadronic energy scale of the calorimeter is ±4.9%. The systematic uncertainty from 
the luminosity calculation is ±1.8%. The error due to the background shape and magnitude 
is found to be by changing the assumption on the vr, K,p composition in the measured 

hadron sample, as described above, as well as by using event samples selected by different 
triggers of the experiment to estimate the background. 

To estimate the theoretically expected cross section, the NLO QCD calculation by Frix- 
ione et al. (FMNR) [§ is used, with rrib = 4.75 GeV and the MRSG [gT]] and GRV-HO ]^ 
structure functions for the proton and photon, respectively. This program supplies the kinemat- 
ics of the b quarks. In order to determine the expectation for the kinematic range where the 
measurement is performed, the b quarks are fragmented to mesons and then allowed to decay 
semileptonically. For the b quark fragmentation the Peterson parameterisation is used with 
£b = 0.006 ± 0.003 [E^. If the fragmentation is done by simply scaling the b quark four vec- 



tor by the factor generated according to the Peterson function the expected cross section in the 
visible range is a{ep bbX — > yuX') = 0.118 nb. However, if the hadronisation of the b 
quark is done together with the generation of a light quark pair then the expected cross section 
in the visible range is a(ep bbX — > yuX') = 0.089 nb. This difference is due to a softer 
transverse momentum distribution of the b hadron in the latter case, which influences the p'^ 
spectrum where a cut is applied. For the theoretical expectation the average of the above values 
is taken and the difference is taken to be the systematic uncertainty. A variation in the renor- 



malisation or factorisation scale, which is defined as fiji = fip = y^l ± pI by a factor of 

2 changes the expected cross section by ±10 %. This uncertainty is added in quadrature to the 
uncertainty due to the fragmentation. The theoretically expected cross section based on a NLO 
QCD calculation is therefore taken to be 0.104 ± 0.017 nb to be compared to the measured 
value of 0.176 ± 0.016 (stat.) to^lj (syst.) nb. The expectation from the FMNR calculation 
in LO for a{ep bbX — > fiX') is 0.069 ±0.008 nb. The corresponding LO QCD expectation 
from the Aroma generator, with which this measurement has been compared previously p5|], 
is 0.038 nb. 

In the Weizsacker- Williams approximation (WWA)p6|] the electroproduction cross section, 
(Tep, is expressed as a convolution of the flux of photons emitted by the electron, /e/^, with the 
photoproduction cross section : 

cTep = a{ep ebbX) = j dyfe/^^a{'jp bbX). 



In order to derive the total photoproduction cross section the measured agp in the visible kine- 
matic range has to be extrapolated to the full phase space. The fraction of the phase space 
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Figure 4: The total photoproduction cross section, cr(7p ^ The horizontal error bar represents 

the range of the measurement. The solid curve shows the expectation of the FMNR NLO QCD calcula- 
tion (full line) with nib = 4.75 GeV and the MRSG and GRV-HO structure functions for the proton and 
the photon, respectively. The shaded area corresponds to the uncertainty if the factorisation scale changes 
by a factor of 2. A change of the renormalisation scale by a factor of 2 leads to a similar result. The 



dashed line represents the prediction of the FMNR NLO QCD calculation if the MRST [ ]280 structure 
function for the proton is used. 



covered by the defined visible range is found to be 12.5%, from the FMNR NLO calcula- 
tion with the hadronization and semileptonic decay, as described above. Correcting for the 
semileptonic branching fractions for a muon originating from a b quark p7| ] and extrapo- 
lating to the full phase space, the total electroproduction cross section for < 1 GeV^ is 
a{ep ebbX) = 7.1 ± 0.6 Ig!? ±1-1 nb. Using the WWA the total photoproduction cross sec- 
tion is a{jp bbX) = 111 ± 10 ± 17 nb averaged over the range 94 < W^p < 266 GeV 
with a mean value of {Wyp) ~ 180 GeV. Here the first and second errors are the experimental 
statistical and systematic uncertainties, respectively. The third is the systematic error due to the 
extrapolation and the branching fraction uncertainties for the semileptonic b quark decays p7| ] . 
The expectation of the FMNR NLO calculation is 63 nb. This measurement together with that 
from EMC and the expectation of the FMNR NLO QCD calculation are shown in figure ^. 



Conclusion 

The open b production cross section has been measured for the first time at HERA using semi- 
muonic decays of the b quarks. The visible cross section a(ep ^ bbX I^X') , in the 
range Q"^ < 1 GeV^ 0.1 < y < 0.8, > 2.0 GeV and 35° < 6^ < 130°, is found to be 
0.176 ± 0.016 {stat.) to^^ (syst.) nb, compared to the expectation 0.104 ± 0.017 nb from 
a NLO QCD calculation. 



From this measurement the total cross section for electroproduction with < 1 GeV^ 
and photoproduction, extrapolated to the full phase space, is calculated to be a{ep ebbX) = 
7.1 ± O.Q{stat.)tli{syst.) nb and a{-fp bbX) = 111 ± 10 (stat.) (syst.) nb at 
an average {W^p) ~ 180 GeV, respectively. The measured cross sections are higher than the 
expectation based on a NLO QCD calculation. 
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